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ABSTRACT 


The effects of wave refraction and damping on swell 
propagation across a wide continental shelf were examined with 
data from a transect of bottom pressure recorders extending 
from the beach to the shelf break near Duck, North Carolina. 
The observations generally show weak variations in swell 
energy across the shelf during benign conditions, in 
qualitative agreement with predictions of a spectral 
refraction model. Although the predicted ray trajectories are 
quite sensitive to the irregular shelf bathymetry, the 
predicted energy variations are surprisingly weak, consistent 
with the observations. The results indicate that small 
amplitude swell is not significantly damped on the shelf. 
However, a large decrease in swell energy levels across the 
shelf (up to 70 %), observed with high-energy incident swell, 
is not predicted by the energy conserving refraction model. 
These energy losses are likely caused by bottom friction. 
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I. INTRODUCTION 


The propagation of swell over a shallow continental 
shelf is a complex process that is still poorly understood. 
While swell can traverse ocean basins with very little loss 
in energy (Snodgrass et al., 1966), strong spatial 
variations in swell energy levels typically occur in shallow 
coastal areas. A variety of processes may affect waves as 
they propagate across the continental shelf from the open 
ocean to the beach. It is well known that refraction and 
shoaling strongly affect the coastal wave climate, in 
particular in areas with complex bathymetry (Munk and 
Arthur, 1951, and O'Reilly; Guza, 1993). Shemdin et al., 
(1980) suggested that bottom friction, percolation through 
the bottom and wave-induced bottom motion, all may cause 
significant damping of swell propagating across a wide 
continental shelf, but the damping rates appear to be 
strongly dependent on the sediment type, bottom 
microtopography and local currents, all of which are often 
unknown. Field measurements show that bottom friction is 
sensitive to the presence of small-scale bedforms (Grant and 
Madsen, 1986), but wave induced ripples are in constant 
transition and extremely difficult to quantify on a natural 
sea bed. 

Very few measurements of the attenuation of swell in 
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shallow coastal waters have been reported. Hasselmann et 
al., (1973) examined the swell attenuation observed during 

the JONSWAP experiment (1968, 1969) with a 160 km long 
cross-shore transect of various instr;iments deployed near 
Sylt, Germany. The observed strong attenuation of swell did 
not agree with generally accepted formulations of bottom 
friction and suggested that bottom damping expressions used 
in wave prediction models are inaccurate or that other 
physical processes such as scattering from small scale 
bottom irregularities cause significant attenuation (Long, 
1973) . 

Young and Gorman (1995) reported similar observations 
from an array of seven instruments spanning the continental 
shelf on the southern coast of Australia. Young and Gorman 
used these measurements in conjunction with the spectral 
wave model WAM [WAMDI Group, 1988] to estimate the 
contribution of bottom friction to the overall observed 
decrease in swell energy across the shelf. Although the 
data analyses span only a period of 18 days at the beginning 
of the experiment (when most instruments were operational), 
and only a few of the instruments were in shallow water. 
Young and Gorman observed significant attenuation of 
energetic swell across the shelf. 

In the present study more extensive observations of the 
transformation of swell over a continental shelf are 
presented. A cross-shelf transect of ten internally 
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recording bottom pressure sensors was deployed offshore of 
Duck, North Carolina extending seaward from the beach to the 
shelf break (100 km from shore), between late July and early 
December 1994. The data set spanned a wide range of 
wind/wave conditions. During periods of light winds, 
remotely generated swell with significant wave heights 
ranging from about 0.1 to 2.5 m were observed. 

Swell spectra observed at each of the instrumented 
sites were compared to predictions of a linear spectral wave 
transformation model initialized with directional wave data 
from a National Data Buoy Center (NDBC), 3-m discus buoy 
located near the seaward end of the transect. A backward 
ray tracing scheme (O'Reilly and Guza, 1993) was applied to 
a high resolution (200 m) bathymetry grid to account for 
refraction effects over the lumpy shelf topography. The 
discrepancies between the energy conserving model 
predictions and measured swell energy levels were used to 
quantify bottom damping effects as a function of the wave 
conditions. 

The experiment, field data and shelf bathymetry are 
described in Chapter II. The effects of shoaling and 
refraction on swell transformation across the shelf are 
illustrated with model simulations in Chapter III. 

Model/data comparisons are presented in Chapter IV, followed 
by a summary and conclusions in Chapter V. 
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II. EXPERIMENT AND DATA ANALYSIS 


The field data used in this study was collected as part 
of the DUCK94 Nearshore Processes Experiment conducted 
offshore of the U.S. Army Corps of Engineers, Coastal 
Engineering Research Center's, Field Research Facility (FRF) 
near Duck, North Carolina. The coast consists of a series 
of relatively straight barrier islands with sandy beaches 
that are fully exposed to the Atlantic Ocean. The 
continental shelf is about 100 km wide and only 20 - 50 m 
deep (Figure 1). The cumulative effect of bottom drag on 
swell traveling across this wide, shallow shelf may cause a 
significant reduction in wave heights on the beaches. A 
transect of wave recorders was deployed on the shelf (Figure 
2) to investigate wave propagation and damping. 

The instrumentation for the experiment consisted of ten 
fully self-contained, battery-powered, internally recording 
bottom pressure sensors deployed along a cross-shelf 
transect extending from the Duck beach to the shelf break 
(Figures 1 and 2, the stations are represented by letters). 
The shallowest instrument X was mounted on a pipe jetted 
into the beach in 6 m depth just outside the surf zone. At 
all other sites (depths ranging from 12 - 87 m. Figure 2) 
the instruments were mounted in the anchor of a surface 
mooring (Figure 3). Heavy steamer chain was used to 
decouple the sensitive pressure sensing instrument from the 
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motion of the surface mooring. The instrument package 
contains a Setra capacitance type pressure transducer, a 
Tattletale microprocessor, and a disk drive for data 
storage. Pressure data was recorded nearly continuously 
with a 2 Hz sample rate during the four-month-long (August - 
November, 1994) deployment. Some malfunctioning data 
acquisition systems were replaced with a cassette tape data 
storage system utilizing a reduced sampling scheme (one 137 
minute long record sampled at 1 Hz every 3 hours). Site B 
suffered significant data loss during the first two months 
of the experiment, and the shallowest instruments X and A 
failed during hurricane Gordon on November 18. The 
shallowest site X (6m depth) and the deepest site I (87 m 
depth) were excluded from the present analysis because the 
beach was not adequately resolved in the numerical 
refraction calculations and high-frequency (> 0.1 Hz) swell 
is strongly attenuated in 87 m depth. 

Measurements of the directional properties of the 
incident swell were available from a National Data Buoy 
Center (NBDC) 3-m discus buoy located within 2 km of site H. 
Although this buoy does not resolve the directional wave 
spectrum in detail, the measurements can be used to 
characterize a mean swell propagation angle and a 
directional spreading factor (O'Reilly et al., 1996). 

Surface height spectra were computed from 12-hour-long 
bottom pressure records using a linear theory depth 
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correction. Relatively long (12 hour) data records were 
used in the analysis because the travel time of the swell 
traversing the shelf is of the order a few hours. Non¬ 
stationary conditions (i.e., temporal variations in spectral 
levels of more than 30% over a 12 hour run) were excluded 
from the analysis because the model predictions do not 
account for time lags in swell arrivals at different sites. 
The analysis was restricted to longer period (0.05 - 0.10 
Hz) waves which are usually remotely generated and feel the 
bottom on the entire shelf. Periods of moderate to strong 
winds (speeds > 10 m/s) with possibly significant generation 
effects at wave frequencies < 0.10 Hz were also excluded 
from the analysis. During the periods of light winds 
considered in this study, currents on the shelf were 
predominantly tidal with speeds generally less than 50 cm/s 
(Haus et al., 1995). The long-period swells considered here 
are not significantly affected by shelf currents. The 
analysis was further restricted to cases with mean swell 
propagation directions (measured near site H) within +/- 35° 
from normal incidence to the shelf break (065° < 0^ < 13 5°) . 
Observations of larger northerly or southerly swell 
incidence angles were excluded because waves approaching the 
shelf at large oblique angles are strongly refracted over 
the continental slope seaward of the instrumented transect, 
and thus the deep water directional properties of these 
waves are not well represented by the NDBC buoy measurements 
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collected at the shelf break (site H). After the various 
rejection criteria were applied, the original data set of 
248 observations was reduced to 71 observations. 

The observed variability in swell energy levels on the 
shelf is summarized in Figure 4 with the total swell 
variances at four sites spanning the shelf. Variation in 
swell energy across the shelf are generally small (< 30%), 
with the exception of a single event from julian days 290 to 
295. During this time frame when maximum incident swell 
energy levels occurred at the shelf break (site H), a 
significant reduction (up to 70%) in energy is observed at 
the shallower sites. 

Accurate predictions of swell refraction requires 
detailed knowledge of the shelf bathymetry. A high 
resolution digital bathymetry database was available from 
the National Ocean Service (NOS), National Geophysical Data 
Center (NGDC). Unfortunately this data base contained large 
gaps extending from 36.2° to 36.8° N, and from the beach 
(75.8° W) to 74.8° W. To fill these gaps, additional 
bathymetric surveys were conducted during instrument 
deployment and recovery cruises with a precision depth 
recorder mounted on the hull of the R/V Cape Hatteras, with 
a track spacing of order 1 km. The fathometer measurements 
were detided using sea level data from a tide gauge located 
near site A on the FRF pier. These corrections are accurate 
only in the vicinity of the tide gauge (i.e., the inner 
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shelf sites A-C) but the errors (< 0.5 m) are negligibly 
small compared to the water depths of the mid- and outer- 
shelf (30-100 m). 

Data from the combined NOS and R/V Cape Hatteras 
surveys was used to create a bathymetry grid for the area 
35° N to 38° N and 74° W to 76° W. The North American Datum 
of 1927 (NAD27) and Mean Low Water (MLW) were used as 
horizontal and vertical references. Surveys with a World 
Geodetic System 1984 (WGS84) horizontal grid reference were 
converted to NAD27 using the Abridged Molodensky Datiam 
Transform Equations (DMA TR8350.2, 1987). A uniform 
bathymetry grid with six second horizontal resolution (-200 
meter) was obtained from the surveys using the Delaunay 
tessellation interpolation method of Watson (1982). This 
method produces a network of near equiangular triangles with 
the vertices being depth soundings and grid points linearly 
interpolated from the plane passing through the vertices. A 
small amount of grid distortion (about 2.5% at the north and 
south ends of the grid) results from transforming a Mercator 
projection graticule to Cartesian coordinates. The maximum 
wave propagation direction errors (»1°) due to this 
distortion are negligibly small (O'Reilly and Guza, 1993) . 

Other potential bathymetric errors that may exist 
include navigational inaccuracies of survey vessels for the 
older surveys and temporal changes in shallow water 
bathymetry due to accretion and erosion of the sandy bottom. 
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These errors are difficult to quantify and may cause 
significant inaccuracies in the shallow water regions where 
small changes in depth strongly affect swell propagation. 

The North Carolina shelf is characterized by a broad 
mid-shelf region with multiple ridge-like features that are 
roughly aligned with the coast line (Figure 1). These 
ridges with amplitudes of the order 5 m (Figure 2), may 
contribute significantly to the refraction of low frequency 
swell propagating across the shelf. 
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III. 


SPECTRAL REFRACTION COMPUTATIONS 


To quantify the importance of the shelf topography in 
the observed spatial variations in swell energy, spectral 
refraction computations were carried out for a wide range of 
deep water incident wave conditions. The incident wave 
field was assumed to be stationary and spatially 
homogeneous, and fully described by a frequency-directional 
spectrum effects of wave generation, 

nonlinear interactions and dissipation on the shelf were 
neglected. For a given deep water spectrum Eo(f,&o) 
predictions of the transformed spectrum Ef'f,0; at eight of 
the instrumented sites A-H (Figure 1) were obtained with a 
backward ray tracing technique described in O'Reilly and 
Guza (1991). 

From each site rays were traced in all possible 
directions back to deep water using the ray equations (Munk 
and Arthur, 1951; LeMehaute and Wang, 1982): 


dX 

dS 


= cos© 


( 1 ) 


dY 


sin© 


( 2 ) 
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dQ 

dS 


1 

C 


sin0 


dC 


cose^ 


(3) 


where C is the phase speed, S is distance along a ray and 0 
indicates the direction of wave propagation. The horizontal 
phase speed gradients were calculated using a second degree 
polynomial fit to the local bathymetry grid, and equations 
1-3 were integrated using a 4*''' order Runge-Kutta method. 

Rays were initially computed for all possible shallow 
water angles 0 at 1° increments and are terminated upon 
reaching deep water, land or the boundaries of the grid. 
These angles were subsequently bisected with additional rays 
until the spacing of the resulting deep water angles 0^ of 
adjacent rays was everywhere less than 2.5° (see O'Reilly, 
1991, and O'Reilly and Guza, 1991, for further details). 

The ray trajectories for a given shallow water site yield an 
estimate of the inverse direction function V : 


0 ^ = r(f,Q) 


(4) 


which defines the deep water incidence angle 0^ as a 
function of the frequency f and the shallow water refracted 
propagation direction 0 (LeMehaute and Wang, 1982). 
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Examples of the inverse direction function F are 
shown in Figures 5 and 6. Figure 5 shows the inverse 
direction functions of sites H (outer-shelf) and C (inner- 
shelf) for a frequency of 0.10 Hz. At site H (49 m depth), 
0.10 Hz swell barely feel the bottom and the F function is 
nearly a 1:1 linear relationship. On the inner shelf (Site 
C), the cumulative effects of refraction over the wide, 
lumpy shelf are evident in the sensitivity of the inverse 
direction function F to the shallow water angle 0. Figure 
6 shows the inverse direction functions of sites H and C for 
0.07 Hz waves. These lower frequency waves sense the bottom 
in deeper water and are already significantly refracted at 
site H. The relatively strong refraction effects noted at 
site C for southerly deep water angles (€>„ > 150°) are 
caused by the curvature of the coast to the south and a 
group of shoals near Cape Hatteras (Figure 1). 

Once the inverse direction function F(f, 0) is 
evaluated for a given site, the frequency-directional 
spectrum follows from the transformation relation 

E (f,9) ^ K 

where K is the wave number, Cg is the group velocity, and 
subscripts indicate deep water values (Longuet-Higgins, 

1957, and LeMehaute and Wang, 1982). Discritizing the deep 
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water spectrum in finite frequency-direction bands 

with widths (Af,A&^), the energy transformation factor 
K(fi,0oj) for each band follows from integrating (5) 


K(f.,e .) = -jrArsr (df fdS ( 6 ) 

i' °j' Af A0 J J K c 

o o g 

where the integration limits include all (f,0) rays that 
terminate within the (fi,0oj) band (i.e., \f - < A.f/2 ; 

irrf,©^ “ < A0^/2) . Finally, predictions of the 

transformed frequency spectrum E(f) at the instrumented 
sites are readily computed for any deep water incident wave 
spectrum by multiplying by the energy 

transformation coefficients K 

E(f.) = E K(f.,e JE rf.,0 J A© 

1 1 oj o' l' OJ o (/) 


To investigate the effects of refraction on typical 
swell arrivals from a single remote source, numerical 
simulations were carried out with a simple cosine power deep 
water directional distribution of energy: 


E f© ; « cos^" 

o ' o' 


© - © 

O MEAN 


( 8 ) 
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The mean propagation direction 9„ean varied from 50° to 
150° (100° is approximately normal to the shelf) and typical 
directional width values of M = 25, 50 and 100 were used in 
the simulations (Figure 7). The swell frequency was varied 
from 0.05 to 0.10 Hz with a finite bandwidth of 0.01 Hz. 

Example model predictions of the transformed swell 
variance at sites A, C, E, and H (Figure 1) relative to deep 
water are shown in Figure 8 as a function of ^ = 

0.05, 0.07 and 0.10 Hz, and M = 50. Although the ray 
trajectories are sensitive to the irregular shelf bathymetry 
(Figures 5 and 6), the predicted swell energy levels are 
generally within a factor of 2 of the deep water value and 
only weakly dependant on Not surprisingly, the 

predicted energy variations on the shelf are very weak for 
relatively short wavelength 0.10 Hz swell (Figure 8c) and 
more pronounced for lower frequency waves (Figure 8a,b). 

It should be noted that very low frequency 0.05 Hz swell 
(unusual at this site) is significantly affected by 
refraction before reaching the shelf break (Site H, Figure 
8a) . For large oblique incidence angles, (<70° and >130°) , 
refraction causes consistent reductions in energy close to 
shore (e.g., sites A and C). 

Small differences (less than 20%) in energy levels 
predicted at site C for different values of the directional 
spreading parameter M (25, 50, 100, Figure 9) indicate that 
the transformation of swell across the shelf is relatively 
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insensitive to the width of the directional spectrum. 
Overall, the model simulations suggest a weak sensitivity of 
swell transformation across the shelf to the deep water 
incident wave conditions. 
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IV. MODEL/DATA COMPARISONS 


Swell spectra on the shelf estimated from 71 twelve- 
hour-long data records (the selection of these records is 
described in Chapter II) are compared here to predictions of 
an energy-conserving spectral refraction model (described in 
Chapter III). The model predictions were initialized with 
directional wave measurements from an NDBC 3-m discus buoy 
located at the shelf break near site H. Estimates of the 
directional distributions of incident swell energy in 0.01 
Hz wide frequency bands (centered at f = 0.05, 0.06, 0.07, 
0.08, 0.09, 0.10 Hz) were extracted from the buoy 
measurements using the Maximum Entropy Method (Lygre and 
Krogstad, 1986). These estimates do not resolve the 
directional spectrum in detail, but, fortunately the 
refraction model predictions are not overly sensitive to the 
directional properties of incident swell (Figures 8,9), and 
mean swell incidence angles are well characterized by the 
buoy measurements (e.g., O’Reilly et al., 1996). Weak 
changes in swell propagation directions between deep water 
and the buoy location (Figures 5,6) were neglected. 

To account for small energy variations between site H 
and deep water owing to shoaling and refraction (note that 
energy levels at site H can differ significantly from deep 
water at low swell frequencies. Figure 8) the swell energy 
spectrum observed at site H was first transformed back to 
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deep water with Equation 6, and subsequently transformed 
across the shelf to all shallower instrumented sites. 

Example comparisons of observed and predicted swell 
spectra and total variances are shown in Figures 10 - 15 
(observed and predicted values are identical at site H, 
where the model predictions were initialized). In the 
majority of the cases analyzed here, with generally benign 
conditions (significant wave heights < 0.5 m), the model 
predictions yield a gradual and weak (less than 25% at the 
shallowest site) decrease in swell energy levels across the 
shelf, in reasonable agreement with the observations (e.g.. 
Figures 10 and 11). The agreement of observed and predicted 
spectra is sometimes poor at low frequencies where energy 
levels are relatively weak (e.g., 0.05 Hz in Figure 10), 
possibly owing to inaccuracies in the buoy measurements or 
the sensitivity of the refraction model to bathymetry 
errors. 

In the relatively few observations with energetic swell 
(significant wave heights > 1.25 m) a large (up to 70% at 
the shallowest site) decrease in spectral levels across the 
shelf is observed but not predicted by the refraction model 
(Figures 12, 13). The observed attenuation is fairly 
uniform over the spectrum (Figure 12) suggesting that the 
effects of nonlinear wave-wave interactions are small. The 
energy losses are likely caused by bottom friction which is 
believed to depend non-linearly on the magnitude of near- 


18 



bottom velocities (e.g., Grant and Madsen, 1986, and Tolman, 
1993) . 

In some cases with relatively small waves, the 
refraction model consistently under-predicts swell energy 
levels at each of the shallower sites (e.g.. Figures 14, 

15). Similar discrepancies were reported by Young and 
Gorman (1995) and attributed to spatial variations in deep 
water incident wave conditions. An alternative explanation 
for these discrepancies is the limited accuracy of the 
incident deep water directional information obtained from 
the NDBC buoy. NDBC buoys are known to over-predict the 
width of the directional spectrum, in particular with low- 
energy swell conditions (O'Reilly et al., 1996), and thus 
cause a bias in the refraction model predictions. 

Comparisons of predicted and observed swell variances 
for all 71 data records are summarized in Figure 16. In 
most cases predicted and observed variances agree within +/- 
30% indicating that the damping of swell across the 
continental shelf (not accounted for by the model 
predictions)is generally weak. However, during the most 
energetic swell arrivals (deep water variances > 10^ cm^), 
the energy levels observed at the shallower sites are 
consistently much lower (up to 70%) than predicted. These 
large discrepancies suggest that large amplitude swell 
propagating across a wide shallow shelf is significantly 
attenuated by bottom friction. 
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V. 


CONCLUSIONS 


The effects of wave refraction (due to depth 
variations) and damping (due to bottom friction) on the 
propagation of swell across a wide irregular continental 
shelf were examined with data collected offshore of Duck, 
North Carolina. A cross-shelf transect of 10 bottom 
pressure recorders, extending from the beach (6 m depth) to 
the shelf break (87 m depth), was deployed for a four-month¬ 
long period spanning a wide range of conditions. 

Measurements collected during periods of strong local winds 
or rapidly changing conditions were discarded to eliminate 
variability owing to local generation and time lag effects. 
The remaining data records selected for analysis generally 
show weak variations in swell energy levels across the shelf 
during benign conditions, but a strong decrease in energy 
from the shelf break to the beach when incident swell energy 
levels were high. 

The effects of the irregular shelf bathymetry on the 
propagation of swell was investigated through simulations 
with a spectral refraction model for a wide range of 
incident wave conditions. Although the predicted ray 
trajectories are quite sensitive to the multiple, ridge-like 
bathymetric features on the shelf, the predicted energy 
variations for realistic swell spectra are surprisingly 
weak. Pronounced refraction effects, evident in a strong 
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decrease in swell energy across the shelf, are predicted 
only for large oblique swell incidence angles. 

Predictions of the energy conserving refraction model 
agree reasonably well with the weak variation in swell 
energy levels observed across the shelf during benign 
conditions. These comparisons indicate that small amplitude 
swell is not strongly affected by bottom friction. However, 
the large decrease in energy levels across the shelf 
observed with high-energy incident swell is not predicted by 
the refraction model. This attenuation is likely caused by 
bottom friction which is believed to depend non-linearly on 
the magnitude of near bottom velocities. 

The model-data comparisons presented here provide only 
a crude estimate of the importance of bottom damping in the 
propagation of swell across the continental shelf. The 
refraction predictions may have significant errors owing to 
the limited resolution and accuracy of the directional buoy 
measurements of incident swell. Detailed measurements of 
incident wave conditions in deep water are needed to obtain 
quantitative estimates of energy losses owing to bottom 
friction and other dissipative processes. 
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APPENDIX 


Figure 1. Bathymetry of the North Carolina shelf. The 
dashed line indicates the instrumented transect. The 
pressure sensor sites are indicated by letters. 

Figure 2. Cross section of the instrumented transect 
(dashed line in Figure 1). The pressure sensor sites are 
indicated by letters. (Site D is 7.5 km south of the 
transect). 

Figure 3. Bottom pressure sensor mooring schematic. 

The instrument package is housed in the anchor and decoupled 
from the surface buoy motions by shock-absorbing steamer- 
chain . 

Figure 4. Observed swell variability on the shelf for 
the 71 data runs analyzed in this study, (a) Swell variance 
at outer shelf site H, (b-d) swell variance at sites E (mid¬ 
shelf) , C (inner-shelf) and A (nearshore), normalized by the 
swell variance at site H. 

Figure 5. Inverse direction function (F) for f = 0.10 
Hz. Upper panel: Site H (outer shelf). Lower panel: Site C 
(inner shelf). 

Figure 6. Inverse direction function (F) for f = 0.07 
Hz. Upper panel: Site H (outer shelf). Lower panel: Site C 
(inner shelf). 

Figure 7. Cosine power directional distribution for 
three values of the directional width parameter M. 


23 



Figure 8. Predicted energy relative to deep water 
versus mean incident wave propagation direction at sites A, 
C, E and H for f = 0.05, 0.07 and 0.10 Hz and a cosine power 
directional distribution with width parameter M = 50. 

Figure 9. Predicted energy relative to deep water 
versus mean incident wave propagation direction at site C 
for f - 0.07 Hz and a cosine power directional distribution 
with width parameter M = 25, 50 and 100. 

Figure 10. Predicted (solid curve) and observed 
(dashed curve) swell spectra at four sites spanning the 
shelf on 12 September, 1994. 

Figure 11. Predicted and observed swell variance 
versus distance from shore, on 12 September, 1994. 

Figure 12. Same as Figure 10 for 18 October, 1994. 

Figure 13. Same as Figure 11 for 18 October, 1994. 

Figure 14. Same as Figure 10 for 5 August, 1994. 

Figure 15. Same as Figure 11 for 5 August, 1994. 

Figure 16. Ratio of observed and predicted swell 
variances at sites E, C and A versus deep water swell 
variance. 
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